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Abstract 


In this article, we formulated a mathematical model for the spread of the COVID-19 disease 
and we introduced quarantined and _ isolated compartments. The next generation matrix 
method was adopted to compute the basic reproduction number (Ro) in order to assess the 
transmission dynamics of the COVID-19 deadly disease. Stability analysis of the disease free 
equilibrium is investigated based on the basic reproduction number and the result shows that it 
is locally and asymptotically stable for Ro less than 1. Numerical calculation of the basic 
reproduction number revealed that Ry < 1 which means that the disease can be eradicated from 
Nigeria. The study shows that isolation, quarantine and other government policies like social 
distancing and lockdown are the best approaches to control the pernicious nature of COVID-19 
pandemic. 
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Introduction 


Mathematical model and its computer simulations are great tools that can be used to 
control and predict infectious diseases globally. The COVID-19 pandemic which is 
popularly known as coronavirus pandemic, a recent pandemic caused by severe acute 
respiratory syndrome coronavirus (SAR-COVID-2). The new virus was named by 
International committee on Taxonomy of viruses (ICTV) as “SARS-COVID-2” because it 
is generally related to the coronavirus responsible for the SARS outbreak of 2003 (WHO 
2019). The recent pandemic was first identified in December 2019 in Wuhan China and 
it was declared to be a Public Health Emergency of International concern. 


On 30 January, 2020 and as pandemic on 11 March 2020 (WHO 2020a). As of 19 
September 2020, the disease was confirmed in more than 30.3 million cases reported 
globally in 188 countries and territories with over 949,000 deaths and over 206 million 
have fully recovered globally (CSSC 2020). 


The pandemic was first confirmed and announced in Nigeria on 27 February, 2020, 
when an Italian citizen in Lagos tested positive for the virus and the second case was 
reported on 9 March, 2020 in Ewekoro, Ogun State when a Nigerian citizen had contact 
with the Italian citizen (NCDC 2020a, Naclean et al., 2020). 


Recently, the spread of covid-19 in Nigeria increases continually as the latest record 
provided by the Nigeria Centre for Disease Control (NCDC) declares Nigeria to reveal a 
total of 56,956 confirmed cases (NCDC 2020b). On 18 September 2020, Nigeria recorded 
221 new confirmed cases and 1 death out of a total of daily test of 2,609 samples carried 
out across the country, the 221 new cases reported from 18 state - Lagos (59), Abia (46), 
FCT (22), Gombe (20), Plateau (17), Rivers (11), Bauchi (7), Benue (6), Ekiti (6) Imo (6), 
Kaduna (4), Kwara (4), Ondo (4), Ogun (3), Osun (3), Bayelsa (1), Edo (1), Kano (1) and 
these results bring Lagos state total confirmed cases to 18,827, followed by Abuja (5,526), 
Oyo (3,226), Plateau (3,192) Edo (2,611), Kaduna (2,326), Rivers (2,220), Delta (1,799), 
Ogun (1,758), Kano (1,734), Ondo (1,594), Enugu (1,234), Ebonyi (1,035), Kwara (1,013), 
Abia (881), Kastina (845), Osu (810), Gombe (799), Borno (741), Bauchi (689), Imo (557), 
Benue (473), Nasarawa (447), Bayelsa (394), Jigawa (322), Ekiti, (313), Akwa Ibom (228), 
Niger (258), Anambra (232), Adamawa (230), Sokoto (161), Taraba (95), Yobe (73), and 
Kogi (3) (Oyekanmi, 2020). 


This work is licensed under a Creative Commons Attribution 4.0 Unported License 
ISSN: 2141 - 4203 Science and Education Development Inst., Nigeria 


Continental J. Biological Sciences 
Agbata et al. (2020) 13 (2): 1 - 16 


The Covid-19 outbreak is seen as the greatest global threat worldwide because it 
spreads globally and millions of confirmed cases, accompanied by thousands deaths 
globally. As of 14 September 2020, WHO reports 28,637,952 confirmed cumulative 
cases with 917,417 deaths in the world and 18 September 2020, total of 56,956 confirmed 
cases, total of 48,305 discharged and total deaths of 1,094 occurred in Nigeria (NCDC 
2020a). 


COVID-19 is spread primarily from human to human through direct or indirect close 
proximity with contaminated surface, objects or with secretion droplets, saliva or 
respiratory secretions expelled from mouth or nose of an infected individual coughs, 
sings, sneezes or speaks. The sign and symptoms of COVID-19 usually start two to 14 
days after contact with infected person or contaminated surface, this period of time 
after exposure and before having symptoms is known as incubation period. Common 
sign and symptoms of COVID-19 include, cough, fever, tiredness, shortness of breath, 
sore throat etc. Older people with underlying medical challenges like serious heart 
diseases (heart failure, coronary arterial disease or cardionryopatty), cancer, type 2 
diabetes, sickle cell disease, high blood pressure (WHO 2020b). 


Mathematical modelling and its computer simulation has attracted the attention of 
many researchers because of its striking applications in the area of prediction and 
control of infectious diseases. Many researchers have proposed mathematical models 
for prediction of infectious disease (Agbata et al., 2019a), Proposed Mathematical 
Assessment of the transmission dynamics of HIV/AIDs with treatment effects. Their 
model comprises of four compartments considering treatment class as the control 
measure. (Tian et al., 2020), proposed a mathematical model for simulating the 
phase-based transmissibility of a novel COVID-19. 


In their study, they developed a Bats-Hosts Reservoir-people transmission network 
model for simulating the potential transmission from source of infection to the human 
infection. The basic reproduction number from reserviour to person analyzed was 
estimated to be 2.30 and 3.58 from person to person. 


Their model revealed that the transmissibility of SARS-COV-2 was higher than the 
middle east respiratory syndrome within the Middle East countries. Our main 
contribution in this article is addition of new compartments such as quarantine, 
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isolation and death which are not considered in most existing models, the addition of 
these compartments give rise to a number of analysis that can be used to _ study 
COVID-19 transmissibility intensively. 


Mathematical Formulation of the Model 


In this section, we developed a mathematical model on COVID-19 based on some 
realistic mathematical assumptions about Covid-19. The total population is subdivided 
into seven compartments, namely susceptible S(t), exposed E(t), quarantine Q(t), 
infected (/,), isolate I,(t), Death D(t) and recovered R(t) let A be the constant 
recruitment rate to the susceptible population, the rate at which people become 
susceptible to the disease and m the rate which the susceptible individuals become 
exposed to the disease, # is the rate of infection. w and a are the rate at which the 
quarantine individuals because susceptible and susceptible individuals are removed 
from the population respectively. The exposed individuals move to the quarantine class 
through the rate b and get recovered through the rate a. 


The isolated individuals get recovered at the rate k and pl is disease induced death 
rate. Quarantined individuals get infected at the rate @ and infected individual 
recovered through the rate d, the natural death is w in every compartment. 


The model is developed under the following mathematical assumptions 
« The recovered individuals do not become susceptible to the disease again. 
* Potion of the infected individuals and isolated individuals can die due to disease. 
« Every susceptible individual has equal chance of being infected. 
« New recruit enters the population through birth and migration. 
* The quarantined individuals can be infected when they come in contact with 
contaminated surface or infected person. 
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Figure 1: Schematic diagram for the model 
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Based on the mathematical assumptions and the schematic diagram in Figure 1, we 
formulate the following differential equations. 


<= hwQ — (m+ a+ py)s 


<= ms — (uy +at+b+ BPE 


dIn 
“t= BE+0Q— (uy tht+dt wn 


<2 = bE- (1, +9 +w)Q 


dl, 

ae Aly — (Hy + Le + k)Is 
A= as +aE + dl, + kl, —'R 
dD 

Gro H2Un + Is) 


N BD OF FP WO NY FR 


Where N(t) = S(t) + E(t) + 1,(t) + Q(t) +: 1,(t) + RW) + D(H) 


Qualitative Analysis of the Model 
In this section we carry out the model analysis such that the equilibrium state, disease 
free equilibrium, and endemic equilibrium in order to determine stability of the model. 


EQUILIBRIUM STATE 

At the equilibrium point we have. 
A+ws—(m+a+t+yp),)s =0 

ms — (u, t+at+b+Pfp)E=0 
BE+0Q—(u,+h+d+p)I, =0 
bE — (wu, +86+@)Q=0 

AL, — (ty + Hz + kD, = 0 
as+aE+dl,+kI,—p~,R=0 
poUn +1.) = 0 


The Disease Free Equilibrium of the Model 
At the disease free equilibrium 


dE dl dl dD dR d 
*=0,@=0,“4=0,7%=0,%=0,.2=0 
dt dt dt dt dt dt 
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The disease free equilibrium of our model equation is 
++ 7p+o+y t+ ptpt)y= A 

(st.Et 1*Qt1,°R*D*) = (—<—*——.0,0,0,0,0,0) 

The Endemic Equilibrium of State 

The endemic state of our model equations is given below 


oo A(Quytatb+B)(u1+8+@) 


~~ whm—[(uy+a+b+B)(uy+O+w)(m+atp4)] 
« — (Wy +0+w)(mt+atps)H*—A 


E 
bw 

iS B(uy+0+w)(m+a+p4)H*bO(m+at+py)H*-2A 
ee wb(uy+h+d+p2) 
Q" _ (mt+atpy)H*—-A 

7 w 
is A+B (uy, +0+@) (m+atpy )H*bO(m+atpy)H*+2Ar 
- bw(uy+h+d+p) 
R= abwH*+a(u4+0+w)(m+atpy,)H* bw(d—k)G* 

7 Hybw 


Where S* = H* and I,” = G* 


The Basic Reproduction Number (Ro) 

In any disease model, the basic reproduction number (Ro) is highly significant since it 
is the most relevant epidemiological distinguishing feature to determine the nature of 
disease. We therefore, use the next generation matrix approach to evaluate basic 
reproduction number of the model. 
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—(u, +at+b+f) 0 0 
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ee 0 0 
(uy thtd+p2) (ui +y2+k) 
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0 0 0 
0 0 0 


The basic reproduction number (Ro) is the spectral radius of the matrix (FV~') which 
is given by 


qcs* 
Ry = —————_ 
N(uith+d+u2) (ui +H2+k) 
z= qca 
0 NiGnta+py—@) (ui th+d+uz) (Wa +H2+h)] 


Local Stability Of Disease Free Equilibrium 
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dV dV dv dV dv dv ada 
no ge ea a a 
Let X =as+aE+dl,+KI,—p,R 

ax ak ak 8 AK AK 
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The Jacobian matrix for the model is given below 


—(m+atpy) 0 0 0 0 00 
m —(u4,+a+b+B) 0 0 0 00 
0 B —G+htd+m) 8 00 
Gp) 0 b 0 —(1+6+) 0 00 
0 0 h 0 —-@itt+k) 00 
a a d 0 k — 0 
0 0 L2 0 29) 00 
Jef -AD = 
—(m+atpy)—-A 0 0 0 0 0 0 
m —(yt+atbt+B)-~A 0 0 0 0 
0 B —(utht+d+tp)-a 8 0 0 0 
0 b 0 —(@+9+a)-2r 0 0 O 
0 0 h 0 —(y t+hg tk)-a 0 O 
a a d 0 k —yw-AD 
0 0 He 0 H2 0 -A 


The characteristic equation is 


{l- (m+ a+py)—-Al+[-Gatatb+p)—-Alt+[-Gate@t+o)-AlJt+[-Gitmt 
=A] + (—w—a] 4-1} 20 


Since all the eigen values are negative, we conclude that the system is locally and 
asymptotically stable around its disease free equilibrium. 
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Numerical Simulation 
We perform numerical simulation of the model in order to compare the results of our 
model with the real data obtained from published articles by WHO and other 
researchers, considering a total of 56,956 confirmed cases in Nigeria. The numerical 
simulation of our model is performed using MATLAB and set of parameter values 
(Table 1). 











Table 1: The numerical simulation of model 
Parameter | Values Description Source 
1. M_ | 0.1818 Rate at which susceptible individuals Agbata et 
become exposed al. (2019a) 
2. B 0.1493 Rate at which exposed individuals Liu et al., 
become infected 2020 
3. H | 0.0752 Isolation rate Agbata et 
al. (2019b) 
4. M2 | 2.5000 Disease induced death rate Estimated 
5. M1 | 0.0400 Natural mortality rate Estimated 
6. b 2.5000 Quarantine rate Estimated 
7. 0.1500 Rate at which quarantined people Estimated 
become infected 
8. A 0.0058 Rate at which exposed people get Ogala et al., 
recovered 2020 
- - 0.0550 Recovering rate of susceptible individuals | Estimated 
, 0.0057 Rate at which quarantined people Estimated 
become susceptible 
11k 0.0780 Rate at which Isolated people get Estimated 
recovered 
122 2.6000 Recruitment rate Estimated 
13.d 0.2000 Recovering rate of infected individuals Estimated 
14.q 1.5000 Per contact rate of infection Estimated 
15.c 0.1000 Probability of getting COVID-19 Estimated 
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qca 


 N[Qnta+py-@) (uy +h+d+pH2) (Mi +H2+k)] 
— 1.5000x0.1000x2.6000 


~~ 20,000[(0.1818+0.055-+0.4—0.0057)(0.4+0.0752+0.2+2.5)(0.4+2.5+0.2+0.00752)(0.4+2.5+0.078)] 
0.3900 


~~ 20,000[0.6311x3.1752x2.978] 
0.3900 


20,000[0.6311x3.1752X2.978] 
0.3900 


~~ 20000x5.9675 
0.3900 


Ro 





119350 
= 0.0000032678 


Graph of Susceptible population against time 
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Figure 2: Graph of susceptible population 


Figure 2, is the graph of susceptible population against time. We observed that the 
susceptible individuals decrease rapidly as time increases showing that few numbers 
will be susceptible to the disease and the infection would be eradicated from Nigeria 
as time progresses. 
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Graph of Exposed population against time 
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Figure 3. Graph of exposed population 


Figure 3, is the graph of exposed population against time, from the graph, the exposed 
population decreases as time increases, the exposed population decreases to zero 
indicating that with time the COVID-19 pandemic can be controlled in Nigeria since the 
exposed population tends to zero. 


Graph of Infected population against time 
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Figure 4. Graph of infected population 
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Figure 4, is the graph of infected population against time, it can be seen from the graph 
that infected population decreases monotonically as time increases which means that 
with time no one will be infected and the disease dies gradually from Nigeria. 


Graph of Quarantined population against time 
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Figure 5. Graph of quarantined population 


Figure 5, is the graph of quarantined population against time, initially quarantined population 
increases rapidly, the initial exponential growth could be as a result of government 
intervention measures such as media cover-age on social distancing and period of lockdown, it 
later decreases due to decrease in susceptible and infected population. 


Graph of Isolated population against time 
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Figure 6. Graph of isolated population 


This work is licensed under a Creative Commons Attribution 4.0 Unported License 
ISSN: 2141 - 4203 Science and Education Development Inst., Nigeria 


13 


Continental J. Biological Sciences 
Agbata et al. (2020) 13 (2): 1 - 16 


Figure 6, is the graph of isolated population against time, it is observed that isolated 
population decreases in Nigeria as time increases this could be as a result of decease in 
infected population and rapid increase of recovered individuals. 


Graph of Recovered population against time 
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Figure 7. Graph of recovered population 


Figure 7, is the graph of recovered population against time, the recovered individuals 
increase greatly as time increases, this could be as a result of effective isolation and 
decrease in rate of infection. 


Graph of Death population against time 
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Figure 8. Graph of death population 
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Figure 8, is graph of death population against time in Nigeria, the population is at 
high rate when the outbreak occurs initially but later decreases with time showing 
that disease induced death tends to zero as time progresses. 


Conclusion 


COVID-19 pandemic is the second world disaster after the World War II. Government 
of Nigeria has taken some policies to control the catastrophic nature of the deadly 
disease in their own capacity. In this article, we have developed and analyzed a 
mathematical model which incorporates isolation and quarantine strategies as control 
measures. Analysis of the model shows that diseases free equilibrium is locally and 
asymptotic stable. The numerical simulation and analysis of the model revealed that 
COVID-19 pandemic can be eradicated from the population. 
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